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Q ! Abstract 

The generalized Chaplygin gas model represents an attempt to unify dark matter and dark 

-C ! 

energy. It is characterized by a fluid with an equation of state p = —A/p a . It can be obtained 

from a generalization of the DBI action for a scalar, tachyonic field. At background level, this 

model gives very good results, but it suffers from many drawbacks at perturbative level. We show 

| that, while for background analysis it is possible to consider any value for a, the perturbative 

\q ■ analysis must be restricted to positive values of a. This restriction can be circumvented if the 

OO ! 

{Sj , origin of the generalized Chaplygin gas is traced back to a self-interacting scalar field, instead of 

the DBI action. But, in doing so, the predictions coming from formation of large scale structures 
reduce the generalized Chaplygin gas model to a kind of quintessence model, and the unification 
scenario is lost, if the scalar field is the canonical one. However, if the unification condition is 

^ • 

imposed from the beginning as a prior, the model may remain competitive. More interesting 
results, concerning the unification program, are obtained if a non-canonical self-interacting scalar 
field, inspired by RastalPs theory of gravity, is imposed. In this case, an agreement with the 
background tests is possible. 
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I. INTRODUCTION 



A lot of effort has been devoted in the last decades to determine the matter-energy 
content of the universe. This effort concentrates not only on the determination of the 
amount of the different components of the cosmic budget, but also on identifying their 
nature. The different components predicted by the standard model of elementary par- 
ticles (baryons, radiation and neutrinos) are believed to be known with high precision. 
Expressed in terms of the ratio to the density necessary to have a flat spatial section (the 
critical density), the recent estimations lead to a baryonic parameter density today given 
by Qbo = 0.0456 ± 0.0016, while the radiation density is given by Q l0 ~ 5 x 10~ 5 [l|. The 
neutrino density depends on the number of neutrino species and their masses, but estima- 
tions lead to values close to that found for radiation. Summing up all these contributions, 
the so-called ordinary matter (including under this concept radiation and neutrinos) give 
a very low mass compared with the critical mass. However, the general features of the 
spectrum of anisotropics in cosmic microwave background radiation imply that the total 
matter density of the universe must be very close to the critical one. In the reference 
using 7-years WMAP measurements, Baryonic Acoustic Oscillations (BA0) and Super- 
nova type la (SN la) data, the curvature parameter is found to be Sl^o = — 0.0057^o;Qogg: 
the spatial section of the universe is essentially flat, and the total density of mass/energy 
existing in the universe must be very close to the critical value. Hence, most of the mat- 
ter/energy of the universe comes from the contributions not predicted by the standard 
model of elementary particles. 

The dynamics of galaxies and clusters of galaxies presents important anomalies, requir- 
ing the presence of a pressureless, non-baryonic component, dubbed dark matter, to be 
explained. The anisotropy spectrum of CMB and the present stage of accelerated universe 
require the presence of another fluid, in principle exhibiting negative pressure, dubbed 
dark energy. For a recent review of these evidences, see reference The main quoted 
candidates to be the constituents of dark matter are axions and neutralinos 3| , but these 
candidates come from fundamental physical theories not yet tested experimentally. For 



dark energy, the most natural candidate is the cosmological constant [4]. An important 
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alternative to the cosmological constant as the dark energy component is a dynamical 
scalar field called quintessence All these proposals to the dark sector of the universe 
have their advantages and their drawbacks, that are exposed in the quoted references. 

A quite different alternative to describe the dark sector of the universe is the so-called 
unified scenario where both components are represented by a unique component playing 
at same time the roles of dark matter and dark energy. The prototype of the unification 
scenario is the Chaplygin gas (CG) model 10] which, in its generalized form, reads 



A 

P= — a> 



where a and A are two free constants. When a = 1 we recover the original Chaplygin 
gas model, which has an interesting connection with the Nambu-Goto action [loj]. Using 
the Friedmann-Lemaitre- Robertson- Walker metric and inserting the equation of state ([1]) 
into the conservation law for a fluid 



p + 3-(p + p) = 0, 



(2) 



(dot means derivative with respect to the cosmic time) we can express p as function of 
the scale factor a: 



P 



A + 



B 



i 

l + a 



(3) 



a 3(l+a) 

where B is an integration constant, indicating that the GCG behaves as pressureless 
matter for a — > and as a cosmological constant when a — > oo. This is the main idea of 
the unification program. 

The Generalized Chaplygin Gas (GCG) model has been tested against many different 
observational data, like type la Supernova(SN la) ll|, Cosmic Microwave Background 
Radiation (CMB) 12Nl4|. matter power spectrum (PS) 15|, |l6|. The general scenario 



emerging from these different tests is not uniform, depending on the priors and on the 
statistics scheme. Initially, it has been argued that the CMB tests [3, [3] favor a scenario 
around a = 0, which essentially reduces the GCG model to the ACDM model. But, in 
reference 14J ]. by using the integrated Sachs- Wolfe effect, models with a ~ 300 seem 



also to be favored. On the other hand, SN la tests favor negative values of a 



The 
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perturbative analysis for negative values of a is not allowed due to resulting negative 
values of the square adiabatic sound velocity [l6| 

c 2 = t = ^- (4) 

We intend to address in the present work the last feature of the GCG model. Is it 
possible to consider negative values of a using a perturbative analysis? Negative square 
sound velocity appears, for a < 0, in the fluid representation of the GCG model. The 
original, more fundamental, representation of the GCG model, is given by DBI action 



17. 



101 ]. which employs a scalar, tachyonic field. In reference [17j it has been shown that, 
for a = 1 (original Chaplygin gas model) the DBI action gives the same perturbative 
expressions as those found in the fluid representation, but that such equivalence, at per- 
turbative level, is broken if a self-interacting scalar field is used to represent the CG model 
instead of the DBI action. 

In the present work, we extend the analysis to the case of the GCG model. It will 
be shown that the complete equivalence of the DBI action and the fluid representation 
remains for a 1. A canonical self- interacting scalar model for the GCG model will 
be developed, and in this case, the equivalence is lost at perturbative level. Hence the 
perturbative analysis can be made even for a < 0. The comparison with the observational 
data will show, however, that positive values of a will be favored, except perhaps when 
the unification scenario is imposed from the beginning. All analysis will be made at 
perturbative level, since the background is the same as in the fluid or DBI representations. 

The canonical self-interacting scalar field, however, hardly can represent a unification 
model for dark matter and dark energy, due to one crucial property: its sound velocity 
is equal to one (in units of velocity of light). Any model for dark matter must have a 
zero effective speed of sound by the end of the radiative era. However, a non-canonical 
self-interacting scalar field, inspired by the Rastall's theory of gravity 18j, may change 
the overall picture. This non-canonical scalar field obeys a dynamical equation very 
different from the Klein-Gordon one, and at same time it can not be reduced to a kind 



of if-essence model, and its effective sound velocity can be set equal to zero [19|. The 



Rastall's theory is based essentially on a modification of the conservation law of the 
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energy-momentum tensor. When, a specific representation for this scalar field is chosen, 
reproducing the background results for the GCG model, the unification scenario leads to 
compatible results even at perturbative level. This is the main result of this paper. 

In the next section, we show the equivalence of the fluid representation with the DBI 
representation for the GCG model at the background level. The equivalence is extend to 
the perturbative level in section III. An equivalent canonical scalar model is developed 
in section IV, and the perturbative analysis is performed in section V. This analysis 
is reproduced for the Rastall's model for the scalar field in section VI, showing that a 
coherent scenario can emerge in this case. In section VII we present our conclusions. 



II. THE DBI FORMULATION 

The DBI Lagrangian is given by the following expression: 

£ = v(7Vi-7yr* (5) 

where T is a scalar (tachyonic) field, and V(T) is the potential for this field. An energy- 
momentum tensor can be constructed from this action. Variation with respect to the 
metric leads to the following expression for the energy-momentum tensor: 

= ~^FlW v = 7l - T iP T» + V ( T Wl-T iP T (6) 

where S is the action constructed from the Lagrangian (jHJ). Comparing with the energy- 
momentum tensor of a perfect fluid, 

Tya, = (ft + p)u^u u - pg^, (7) 

and using the homogeneous, isotropic FLRW metric 

ds 2 = dt 2 - a(t) 2 ^ ij dx i dx j , (8) 

we obtain the following expressions for the density and pressure: 

V(T) 



Pt= t^== , Pt = -V(T)V1-T*. (9) 
Vl-T 2 
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This leads to the equation of state 



Pt 



V{Tf 
Pt 



(10) 



For V(T) = y A = constant, the equation of state (TTOT) represents the traditional Chap- 
lygin gas model. When V(T) is not a constant new possibilities are opened, and the 
richness of the model has been pointed out in reference 201 ] . 
The GCG, characterized by the equation of state 



A 

P = 7 



can be obtained from the action 

C T = V(T) 



(d p Td p T) 



l + a 
2a 



;i2) 



where V(T) = = constant. In fact, from ( fl2l) . we have the energy-momentum tensor 



= VjTHd^Ti^d^T + V(T) 



1 - (d p TdPT) 



l + a 
2a 



(d p Td p T)^ 



9 flu- 



(13) 



A direct comparison with the perfect fluid energy- momentum tensor (I7j), leads to 

V(T) 



Pt 



{d p TdPT) 



l + a 

2a 



1 5 
l + a 



PT = —V(T) 



1 - (d p Td p T)^ 



u, 



yJdpTdPT' 



(14) 

(15) 
(16) 



III. PERTURBATIONS OF THE GENERALIZED DBI ACTION 

The perturbation of the perfect fluid energy-momentum tensor ([7]) is now carried out 
using the synchronous coordinate condition = 0. The components of the perturbed 
energy-momentum tensor are 

ST 00 = Sp, (17) 
ST 0i = {p + p)5ui, (18) 



5T 7 



ij 



-phij - 8pg. 



(19) 



The perturbation of the generalized DBI action ( 112]) leads, on the other hand, to the 
following expressions when V(T) = V = constant: 

1 vf^sf 



5T ( 



oo 



a 



2+a 1 

1 + Q, 



sn 



VT°di6T 



1 5 
1+ct 



VT«5T TA/ _ X 

— i-<to + v(T) 



1-T J 



1 -T~ 



(20) 



(21) 



(22) 



From the general form of the density, pressure and velocity obtained from the energy- 
momentum tensor of the DBI action, (|14|I15|I16)) . we obtain, in the synchronous coordinate 
condition, 



5 p 7 



1 V(T)fa6f 



a 



T 



2 + a 1 
1 + a 



V(T)f«6T 



1 5 
1 + a 



(23) 



(24) 



(25) 



Inserting the expressions (j23|24)25p into (j20|21|22p . we obtain ( |17|18|19| . Hence, the 
perturbation of the generalized DBI energy-momentum tensor is perfectly equivalent to its 
fluid counterpart. This equivalence has been verified for the first time (to our knowledge) 
in reference |2l|]. It is direct to verify that, 

6p T 



6 pa 



Pt 

Pt' 



(26) 



as far as the potential is constant. The case were V(T) is not constant has been analyzed 



for some specific forms of the potential, in reference [21]. In this sense, a negative implies 
a negative squared sound velocity. This can be verified taking the ratio between (124"P and 
( |23l) . In this sense, the extension of the perturbative analysis for negative values of a is 



forbidden, in opposition to what happens with the background tests, like SNIa or BAO, 
where the problems connected with the sound velocity do not appear. 



IV. SCALAR MODEL 

In order to exploit the possibility of a < we must abandon the DBI framework. The 
most general framework besides the DBI is the self-interacting scalar field. If this self- 
interacting scalar field is implemented in the canonical way, there is the mentioned problem 
with the sound speed: cj. = 1. Moreover, while the DBI action is complete equivalent 
to the Chaplygin gas model, the self-interacting scalar field presents many challenges, as 
non-unicity and stability of trajectories, requiring specific initial conditions jzj). For the 
moment, we ignore all these problems and consider the canonical self-interacting scalar 
field as an effective model which asks for a fundamental description before the decoupling 
between matter and radiation, as well as particular initial conditions. This possibility 
has been touched in \x^, but fixing a — 1. But, at least for this case, it has been found 
that neither the unification scenario H] nor the anti- unification scenario have been 



resulted from this self-interacting scalar field approach: the matter density parameter was 
in fact very close to that predicted by the ACDM or the quintessential model. The aim 
of the present investigation is to verify the predictions when all possible values of a are 
considered. 

The action is now given by 

1 



-g 



R - <p, p fi p + 2V{<f>) 



+L m , (27) 



16nG 

where V(4>) defines the potential and L m is the matter Lagrangian. The field equations 



arc 



Rpv - \qiwR = BirGT^ + (p.^-u - ^gf,u4>; P 4>' p + 9p,V (</)), (28) 

= 0, (29) 

□0 = -V4,, (30) 

where the subscript <fi indicates derivative with respect to the scalar field. 
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( a 


3| 






\a 




( a 






a 


V a 



The equations of motion are 

2 ; 2 
= 8vrGp = ^- + ^(0), (31) 

2 ;2 
= -87rG^ = -|- + y(0). (32) 

In the generalized Chaplygin gas model, the energy density behaves as 

Pc = Pc g{a)^, g(a) = A + ■ (33) 



This expression may be rewritten as 



Q c = n c0 g(a)^, Q c = orr2 c , f2 c0 = 2 C ■ (34) 



SnGpc 8nGp c0 

Making the substitution 



3# 2 ' LU 3H* 



H 



we may write 



Hence, we obtain 



^+V(^) = 3£W^, (36) 
^ - V{cf>) = -3flco^-^. (37) 



f i - 1 1/2 
= v 7 ^^^ - Ag-Th I , (38) 

Vtf) = ^.Sgjh+Ag-^X. (39) 



This scalar field model reproduces exactly the background of the generalized Chaplygin 
gas in presence of baryonic matter. When baryons are absent, we find the potential of 
reference Q. 
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V. PERTURB ATIVE ANALYSIS OF THE SCALAR MODEL 



In order to perform a perturbative study of this model, we rewrite the Einstein's 
equations in the following form: 

R»u = 8tvg(t^ - ^g^Tj+cf).^ - g^V{<f>), (40) 
Dcf> = -V^), (41) 
T^ ;M = 0. (42) 

Again, we choose to employ the synchronous coordinate condition h^Q = in carrying out 
the perturbative study. Since all relevant modes are well inside the cosmological horizon 
the final results do not depend on this choice. 

The standard perturbative calculation using the synchronous coordinate condition 
leads to the final set of equations: 

Q 

5 + 2-5- -Q m 5 = 2<p\ - 2V (t> X, (43) 

(X — 

' A = (f)5. (44) 



A + 3-A+ 
a 



a 2 



+ V 6 , 



In these equations, 5 = 5p/p is the density contrast, A = 5<fi, the parameter k is the 
wavenumber resulting from the Fourier decomposition of the perturbed quantities, and 
fl m = Qmo/a 3 is the density parameter for the matter component, Q m o being its value 
today. This density parameter contains two terms, the baryonic component and the dark 
matter component Q m0 = Q^o + ^dmo- Equations ( I4"3"f4"4"|) refer to the Fourier transform 
of the fundamental quantities, that is, we should more properly write 5^ and A&. 

In order to carry out the comparison with observations, it is more convenient to re- 
express equations (I43f44l) using the scale factor as dynamical variable. Moreover, we 
divide both equations by Hq, the Hubble parameter today. The final (dimensionless) 
equations are the following: 



a f 



(klo\+ V 
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A = y, (46) 



where the primes mean derivative with respect to a and Iq = 3.000/i 1 Mpc is the Hubble 
radius. Moreover, the following definitions were used: 



f(a) = ^j^f + n c {a)a\ (47) 

0(a) = v / 3^V / ^ (a)1/(1+Q) " Ma)- a/(1+a) , (48) 
V(a) = hl c0 (gia) 1 ^ + i<7(a)~ a/(1+a) ) , (49) 

VM = ^V'(a), (50) 

V»lfl) = ^Vj(a), (51) 
n c (a) = n c0 g(a)^ 1+a \ (52) 

In these expressions, f2 c0 is the density parameter for the generalized Chaplygin gas model, 
which obeys the flat condition Q c0 + Q m0 = 1 and g(a) is given by ( 1331) . 

Now, we compare the model with the power spectrum observational data from 
the 2dFGRS compilation in the range that corresponds to the linear regime, that is 
0.01 Mpc" 1 < kh" 1 < 0.185 Mpc" 1 . This compilation amounts to 39 data. In the numer- 
ical computation, one important aspect is how to introduce the initial conditions. We use 



the BBKS transfer function 23], which gives the spectrum today for the ACDM model, 



integrating it back to the redshift z = 1000, where the initial condition are fixed. The 



general procedure is described in reference 24]. In comparison with this reference, there 



is an important modification. When a < — 1, the general behavior is the dominance of 
an accelerated phase in the past, approaching a dust dominated universe in recent times. 
This scenario may bring problems concerning the formation of structure (to which extent 
it will be verified in what follows), but may be interesting to implement the idea that the 
acceleration of the universe is a transient phase. In fact, in reference j^, it has been 
pointed out that the acceleration of the universe is slowing down. The important point is 
that, when a ~ —3, the resulting scenario is very similar to the ACDM model until very 
recent times: for a ~ —3 it is found <fi ~ and V(<p) = constant until z ~ 20. The general 
behavior of <fi in terms of z and a is illustrated, for some specific cases, in figure (CQ). For 
a > — 1 this general behavior is inverted. The fact that for a negative enough the model 
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FIG. 1: Behavior of <p as function of a and z, fixing f2 m o = 0.3 for A — 0.1 (left), 0.5 (center) and 0.9 
(right). The darker colors indicate an almost zero value for <j>. 



is essentially ACDM until very recently does not allow to imposed the initial conditions 
using the transfer function with Q m o ~ 0.3 and Qao = 0.7 as it has been done in 24]. 
For a ~ —3 the initial conditions must be imposed, due to computational reasons, about 
z = 20. Because of that, we must take into account the fact that the transfer function 
depends explicitly on the mass parameters. 



The fundamental quantity to be evaluated is the matter power spectrum 

Pk = Si (53) 

which is the Fourier transform of the two-points correlation function of matter distribution 
in the universe. We perform a Bayesian statistical analysis, using first the \ 2 function 

n ( po _ pt \2 

x 2 = E — 2 ' - ( 54 ) 

1=1 0i 

where P%. is the i th power spectrum observational data, with of observational error bar, 
its corresponding theoretical prediction. From this quantity, we construct the Probability 
Distribution Function (PDF), 

P = Ae~ x2/ \ (55) 

where A is a normalization constant. The PDF depends in general, for the GCG model, 
of 5 parameters: A, a, Qdmo ^ c o an d Qko, where the last quantity is associated with 
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the curvature of the spatial section. The power spectrum is expressed as function of h, 
reducing the number of parameters to four. Imposing the flat condition Qko = 0, implies 
only 3 free parameters. We will perform a two-dimensional computation, fixing A. First, 
the computations will be made considering three values of this quantity: A = 0.1, 0.5 and 
0.9. This will allow us to evaluate the influence of this parameter in the final results. We 
will plot the two dimensional PDF as well as the corresponding one dimensional PDF, by 
marginalizing (integrating) over the one of the variables. 

The results indicate a very general pattern. Leaving a free, we remark the appearance 
of the two plateau in the PDF: one for a ~ —2.5, and the other for a ~ 1. The plateau 
corresponding to the positive values of a is higher than the plateau corresponding to 
negative values. Hence, the model predicts, in some sense, that a must be positive. From 
figures ( l2"f3"PJ) it is clear that the PDF is concentrated in the region of positive values for a 
and values for Qdmo in the range 0.2 < fldmo < 0.3. These regions are displaced for larger 
positive values of a as A increases. There is another region with smaller probabilities for a 
negative - see figures (T2T3H1) - which is almost not affected by changing A. These behaviors 
are confirmed by considering the one dimensional PDF for a displayed in figure ([3]). 
Remark that there is a peak near a = — 1 which corresponds to a linear, constant, relation 
between pressure and density. We may attribute this feature to effects corresponding to 
the numerical computation near the singular point - however, even increasing strongly 
the precision of the computation (implying increasing the computational time), this peak 
does not disappear. 

For the dark matter density parameter Qdmo, the PDF is peaked quite generally around 
Qdmo = 0.23. This is a result very similar to the ACDM model, and all framework is very 
close to the quintessence model (whose precise prediction depends on the choice of the 
potential for the scalar field). The parameter estimations at, for example, la, 2cr, etc. for 
the parameter a becomes very doubtful due to the existence of the two plateaus. It could 
be done for the parameter Qdmo due to its almost gaussian probability distribution, but 
even a visual analysis of figures (jlj) shows that they change very little with the increasing 
value of A, remaining always peaked around 0.23. 
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FIG. 2: Two dimensional PDF for the parameters ridmo and a for A — 0.1 (left), 0.5 (center) and 0.9 
right. 




It can be expected that the results may change by imposing some prior. One possibility 
is to restrict a > — 1. This restriction does not change the results, which are displayed in 
figure (pj]). Again, the maximum of probability for Qdmo is sharply peaked around 0.23, 
while there is a plateau for the PDF of a which, in this case, by fixing A = 0.5, begins 
at a ~ 8. However, some interesting effects appear if we fix, from the beginning, the 
dark matter component equal to zero. This amounts to impose the unification scenario 
where the dark matter and dark energy are both represented by the Chaplygin gas model. 
In this case, two free parameters are considered: A and a. In figure ([6]) the results are 
displayed. The main feature to be remarked is that, now the maximum for PDF occurs for 
negative a. There are strong oscillations near a = —1. In this case, the minimum for y 2 
may depend on the precision of the numerical evaluation. It can vary from, for example, 

X 2 ~ 0.30 for a 0.85, to x 2 = 0.34 for a = -0.88 (using SN la, Gold sample, the 

favored value is a = —0.10 [111]). When the unification scenario is not imposed from the 
beginning (the previous cases), it is found typically Xmin = 0.30, in a more stable way. For 
a comparison, for the ACDM, we have Xmin ~ 0.38, but with just one free parameter when 
the spatial section is flat. If the Akaike Information Criteria, AIC, is used, which allows 
to compare models with different number of free parameters (AIC = 2 x N + x 2 taii 
where N is the number of free parameters), the ACDM model remains the best one, 
but the difference is very small: AIC ~ 18 for the scalar GCG model with three free 
parameters, AIC ~ 17— 18, for the unified scalar GCG model with two free parameters, 
while AIC ~ 17 for the ACDM model. 



VI. THE RASTALL'S UNIFICATION SCALAR MODEL 

The results described above are interesting but suffer from a major drawback: a canon- 
ical self-interacting scalar model hardly can represent dark matter, since its sound speed 
is different from zero, in fact it is equal to 1. One of the main properties of dark matter is 
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FIG. 5: Two and one dimensional PDFs restricting a > — 1 and fixing A = 0.5. 





FIG. 6: Two and one dimensional PDFs restricting QdmO = 0. 



a zero sound speed already in the radiation era (at least in the end of the radiation era). 
Hence, the previous model has to be seen as an effective model valid during the matter 
dominated era, which must be complemented for previous era. However, there is another 
possibility to describe a dark matter scalar model, or even a dark matter/dark energy 
scalar model, which is a particular case of the scalar-tensor description of the Rastall's 
theory of gravity. This possibility is explored in some detail in reference 19]. Let us 
follows this approach now. 

Rastall's theory has been proposed in the beginning of the 70's in order to take into 
account a fundamental idea: the usual conservation law for the energy-momentum tensor 
has been tested only in flat space-time, and in curved space-time its employment is a 
pure extrapolation. Hence, Rastall has proposed that the Einstein's equations should be 
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written as 

R P u - -zQ^R = 8nGT^, (56) 

where A is a free parameter. If A = 1, General Relativity is recovered. The Rastall's 
equations may be re-written as, 



Rfa> ~ ifiiivR = 87rG 



2-7 



(57) 



T- ; „=X^V, 7 =|^ (58) 

One of the main theoretical drawbacks of the Rastall's theory is the absence of a La- 
grangian formulation. However, a possible action principle has been proposed in reference 
26j , but outside the framework of riemannian geometry. In reference [23] Rastall's theory 
has been tested against the 2dFGRS data for power spectrum. In order to do so, it has 
been considered a two-fluid model, one obeying the usual conservation law, and the other 
the Rastall's conservation law. This hybrid model was necessary in order to represent the 
baryon component, which must have zero pressure, clustering to form the observed struc- 
tures. The second component should represent dark energy, and can have a zero pressure 
but behaving effectively as a component with negative pressure. The hydrodynamical 
representation leads to a very stringent limit on the parameter 7, reducing it effectively 
to the Einstein's theory. 

Hence, a scalar model was conceived in order to represent the " Rastall's fluid" . In this 
case, the equations read now, 

R»u - \g^R = SnGT^ + <P,^, U - l^Lg^tf,.^ + g^(3 - 2j)V((f>), (59) 

v P vy + (3 - 2 7 )n = (i - i) r ftr - ( 6 °) 

As before, represents the baryons and obeys the usual conservation law. Again, the 
problem of a action principle leading to equations (I59|60p must be considered, and it can 
be treated in a similar way as the formulation of reference 26(. It is worth to mention 
that the modified Klein-Gordon equation fl60l) is similar to some terms appearing in the 
Galileon's theory 28|, but it seems that no exact correspondence can be established 



between the scalar Rastall model and the Galileon framework. 
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The interesting point about this model is that the Rastall's self-interacting scalar field 
can have a zero pressure behavior at perturbative level when 7 = 2. This has been point 
out in reference 19j, without an explicit mention to the Rastall's theory. 

Hence, to have a zero sound velocity, the scalar model must obey the following equa- 
tions: 

Rfiv - ^g^uR = 8nGT pu + <p^(j)- u + g^ u V(4>), (61) 

V p + V <j> + V f ™° = 0, (62) 
<p. a (p' 

where, just for future convenience, we have made the redefinition V(4>) — » —V{4>). In this 
case, the energy-momentum tensor for the scalar field is given by, 

T+, = </>„</>.,* + g liV V{<f>). (63) 

Hence, the energy and the density of the scalar field is given by the following expressions: 

p = 2 + y(0), (64) 

= ~V{(f>). (65) 

Now, we must to proceed as in the previous section to fix <f> and V(<f>). Imposing that the 
scalar field must represent the GCG, we find 

0(a) = v / 3?^ V /(7(a) 1 /(i+«) - Ag(a)-^+ a ), (66) 
V(a) = 3n c0 Ag(a)- a ^ 1+a \ (67) 

The perturbed equations read now, 



5 + 2-5- ~-±6 = # - Vfl, (6* 



a- 3Q 
'—5 

a 2 a 3 

2f + 3-*+(^+y^W = 4>5, (69) 
a [a z J 

where \I/ = 5(f) and 5 is the density contrast of the matter component, as before. Using 
now the scale factor as the variable, the above system of equations take the following 
form: 

rJl + mh = (70) 

[a f(a) J 2a 3 f 2 (a) f 2 (a) 

^i+'SM^m}*-*'- (71) 
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The function / has the same form as in the preceding section, but now <p an d V{4>) are 
given by ( I68|lt39"j) . 

The model is tested against the 2dFGRS power spectrum data. We implement the 
initial condition as described in the previous section. The same considerations can be 
made, grosso modo, about the behavior of the potential of the scalar field. The statistical 
analysis is made in the same lines as before. First, we leave two free parameters, a and 
fldmo and fix A. When A = 0.1, 0.5 and 0.9 the results are displayed in figures (171 IH1I91 . For 
a the results are similar to those found previously with the canonical scalar field. But, for 
fldmo, there are now two maximum of probabilities, one when the dark matter component 
is absent, and the other when dark matter is the only component besides baryons. The 
second maximum is higher than the first one. In all these cases, the AlC parameter is 
about 20. If, as before, a is restricted to be less than —1, the results are essentially the 
same, including for the AlC parameter. 

If the unification scenario is imposed from the beginning, leaving A and a as the free 
parameters, the behavior for a remains the same, while the most probable value for A is 
near A = 0.004 - this leads to a an almost CDM scenario, with a cosmological term very 
near zero. This result can modified as the numerical precision is increased but A = 
(exactly CDM model) seems not to be the favored case. Again, the AlC parameter is 
about 20. 

With respect to the ACDM model, the difference in the AlC parameter is AAIC ~ 3. 
Following the Jeffrey's scale this give a moderate support in favor or the ACDM compared 
with the Rastall's scalar model for the generalized C haplyg in gas {29), 0]- Perhaps, due to 



29. 



30j, it would be more appropriate 



details concerning the implementation of this scale 
to state that both models remain competitive as far as the structure formation test is 
concerned. 
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FIG. 7: Two dimensional PDFs for the Rast all's scalar model with A = 0.1, A = 0.5 and A = 0.9. 




FIG. 8: One dimensional PDFs for the parameter a in the Rastall's scalar model with A — 0.1; A = 0.5 
and A = 0.9. 




FIG. 9: One dimensional PDFs for the parameter tldmo in the Rastall's scalar model with A = 0.1; 
A = 0.5 and A = 0.9. 
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FIG. 10: Two and one dimensional PDFs restricting a > — 1 and fixing A — 0.5 for the Rastall's scalar 
field model. 




FIG. 11: Two and one dimensional PDFs restricting £ldmO = for the Rastall scalar field model. 



VII. CONCLUSIONS 

We have investigated here the possibility that the Generalized Chaplygin Gas (GCG) 
model may be represented by a self-interacting scalar field, instead of using the fluid rep- 
resentation or the Dirac-Born-Infeld (DBI) action for a tachyonic fluid. The main reason 
for this investigation is that both the fluid or DBI representations forbid to extend the 
analysis of the CGG model to negative values of a at perturbative level, due to a negative 
squared sound velocity, which drives strong instabilities. This represents an important re- 
striction, mainly when it is taken into account that some background observational tests, 
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like SN la, favor negative values of a. 

We have shown initially that, for the GCG model, the fluid representation is equiv- 
alent to the tachyonic representation of the generalized DBI action, at background and 
perturbative levels, extending the results of reference Hence, in order to consider the 
possibility of negative values for a in a perturbative analysis, we must use another frame- 
work, which we have chosen, first, that one represented by a canonical self-interacting 
scalar field, and later by a non-canonical scalar field inspired in the Rastall's theory of 
gravity. In these cases, negative values for a do not lead to any problem connected 
with the sign of the (squared) sound velocity. When baryons are presented, there is no 
closed form for the potential representing the GCG in both scalar models, but an implicit 
expression, using the scale factor as variable, can be obtained. 

We have made a statistical analysis using the perturbative expressions for the self- 
interacting scalar field model for the GCG, computing the matter power spectrum for the 
matter component, and confronting it with the 2dFGRS observational data. In principle, 
we allow the presence of dark matter besides baryons and the scalar field. For a canonical 
scalar field, the constraints obtained indicate that a kind of quintessence scenario emerges, 
with Qdmo ~ 0.23 (like in the ACDM model). High positive values of a are favored. 
The best fitting configuration implies typically Xmin ~ 0.30, better than in the case of 
the ACDM model, for which Xmin ~ 0.38. An interesting scenario emerges when the 
unification scenario is imposed from the beginning, fixing Qdmo — 0. The best fitting is 
Xmin ~ 0.30, still smaller than in the ACDM case, but now for a negative. 

The above described scalar model has an important drawback: the velocity of light 
is equal to one, and it cannot represent dark matter at the end of the radiative era, 
and it must be seen as an effective model valid during the matter dominated era, asking 
for a coherent description for the previous era. Hence, another model was introduced 
based on the scalar field version of the Rastall's theory. For this case, the velocity of 
the sound can be set equal to zero, independently of the potential term chosen. For 
this non-canonical scalar field, the quintessence model is not anymore favored and, when 
the dark matter density parameter is free, there are two maximums, one for a universe 
entirely filled by dark matter (in agreement with the power spectrum results for the GCG 
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in the DBI representation), and the other when dark matter is absent (in agreement 
with SN la results). The best fit is achieved, generally, with xLin ~ 0.38, similar to the 
KCDM case, but with one more free parameter. Due to this, and employing the AIC 
criteria and the Jeffrey's scale to evaluate the evidence in favor of one model or other, 
supports marginally the AC DM model, but without a decisive result: the scalar field 
model remains competitive. All these results are valid even when the unification dark 
matter/dark energy prior is imposed from the beginning. 

For the parameter a there is a qualitatively agreement with the SN la analysis. Perhaps 
this may indicate a new concordance model, different from the ACDM. However, a more 
deep perturbative analysis must be performed, mainly using the Integrated Sachs- Wolfe 



effect (a very delicate test for the unified models [3l|, |32|), or even the full anisotropy 



spectrum of CMB. We hope to present this analysis in the future. 
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